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Abstract 

We consider an alternate fourth-order gravity Lagrangian which is 
non-analytic in the Ricci scalar, and apply it to a Robertson- Walker 
, metric. We find vacuum solutions which undergo power-law Inflation. 

Once matter is introduced the theory behaves very much like ordinary 
General Relativity; except that the radiation evolution a ~ i 1 / 2 is not 
OS ■ allowed since it corresponds to R = 0. We comment on the possibility 

*q ' of wormhole solutions in such a theory. 
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1 Introduction 



Inflationary universe models can be constructed by adding higher order cor- 
rections to the usual Einstein-Hilbert Lagrangian for gravity. The first such 
model was constructed using the trace anomaly of conformally coupled fields 
which is due to a term ~ R 2 logR/ //(// = renormalization scale) in the La- 
grangian[l]. However, typical numbers and types of field in the early universe 
would not produce enough inflation, and any inflationary phase would finish 
too abruptly causing an inhomogeneous universe. 

The model was modified by adding a term eR 2 to the Lagrangian which 
results in a quasi-de Sitter expansion [2-4]. The Hubble parameter decreases 
slowly for large e before going into an oscillation phase which can reheat the 
universe [4, 5]. The general Lagrangian so far considered has the form f\ 

C = R + eR 2 + (3R^ v (1) 

The coefficients e and j3 have dimensions of [length] 2 and (in order for stability 
or to prevent tachyonic solutions) certain signs have to be taken. The (3 term 
is only relevant for anisotropic models since in the isotropic case it just alters 
the e coefficient. New gravitational waves have also been found in such 
Lagrangians[6]. 

Various other Lagrangians such as R + eR 2 + f3R? and exp(XR) have been 
considered and generally give power law inflation [7-9] : instead of exponential 
expansion there is typically an expansion law a ~ t m . So long as m > 1 and 
expansion occurs for sufficient time the usual problems, which are solved by 
inflation, can still be solved. 

Some papers have appeared concerning the stability and solutions of gen- 
eral Lagrangians f{R) [11-14]. It is sometimes claimed that such theories or 
their quasi-de Sitter solutions are unstable [11,12] but at least classically this 
is not the case. This was analyzed in detail in ref. [13] for the R + eR 2 theory 
and also found to match ordinary gravity as e — > 0. The so-called small 
perturbations considered in ref. [11] are actually abnormal: they correspond 
to a large dH/dt (if=Hubble parameter). In analogy with a scalar field in- 
flationary model this H term is like a velocity causing you to shoot up the 
scalar potential - but it is eventually damped [15]. 

3 Some works use the equivalent in 4 dimensions R + aR 2 + bC 2 Lagrangian, where C 
is the Weyl tensor. 
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A recent problem has arisen in the case of Bianchi type-I models for the 
Lagrangian (1) for (3 7^ so including the sign that is stable in the absence 
of anisotropy. The model is found to collapse due to the anisotropy [16]: 
this was also found by one of us previously using the equations presented in 
ref.[17] but not taken too seriously at the time since it conflicted with ref. [18]. 
Note that this is more restrictive than for the case of scalar field inflation, 
where spatial gradients can cause collapse[19] but not simply anisotropy [20]. 

Instead in this paper we consider an alternative higher order derivative 
theory first proposed in ref. [21] 



Because the parameters are now dimensionless such a Lagrangian could 
help incorporate higher derivative theories into variable G ones (eg. Brans- 
Dicke[22]) or those involving phase transitions c.f. induced gravity [23]. 
These terms in the Lagrangian might result when quantum gravity correc- 
tions are made. Although we are not aware of them being invoked from string 
theory, which is usually taken to give corrections of the form Lagrangian (1), 
this might change when the theory is better developed or is superceded. 

In any case such a non-analytic Lagrangian is interesting in its own right. 
It is considerably more complex than Lagrangian (1) and requires us to look 
for solutions with R 7^ 0. If it is found to have vacuum Inflationary solu- 
tions, it would be a further extension to the generality of higher derivative 
theories in having such Inflationary behaviour. Since we would expect that 
such corrections to the Lagrangian would be present in the early universe 
,when quantum gravity effects were dominating, they could help explain any 
inflationary epoch. 



C = R + XR^RT/R + t R aM& R a ^ s /R. 



(2) 



2 Field Equations 



We consider an action of the form 




The field equations follow from variations of this action leading to 



= G% + AGT + tG% 



8ttG 



(4) 
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where 



QHV = RVV _ i gilv ^ 



(5) 



= g ^ {R^/R), aT + {BT/R)./ - <T (i2a/3 R aP /R% T 

- [{R^/R).% + {IT/R).%\ + (R af3 R aP /R 2 )r 

+ (i? a/3 iT^/i? 2 ) - 2/T 7 i^/i? + | R a p R a(3 /R, (6) 

Gf = A{R^/R). aT -g^ (R afSl5 R a ^ 5 /R% T 

+ (R aPl5 R^ S /R 2 Y U - \ R^r P R urjrp I R 

+ BT R aPlS R a ^ s /R 2 + \ R af3y5 R a ^ s /R . (7) 



Since R^ v and R a ^ d include second order derivatives of the metric, the field 
equations are fourth-order partial differential equations, which is the same 
order as the normal higher order gravitational theory. 

Before considering a Robertson- Walker metric we can find the equation 
for a maximally symmetric spacetime. Such a spacetime has 



Substituting (8) and (9) into (6) and (7), respectively, and since the covariant 
derivatives of g^ v are equal to zero, we have 




(8) 



Raf3y5 = —R (<?a 7 9/38 ~ 9a8 #7/3) /12. 



(9) 



Gf 



= BT R af3 R aP /R 2 - 2 RT R V JR 
+ \trR aP K*P/R=±g'»'R, 



(10) 



= -2 R\ rp R vaTp I R + FT R aPlS R a ^ 5 /R 2 
+ \ <T Raw R a ^ s /R =±R <T . 



(11) 



Substituting (10) and (11) into (4), we find 



\ g" u R(l - A/4 - r/6) = 8nG/c 



(12) 



i2=-5-T/(l-A/4-r/6) 



(13) 
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which is just the Einstein's field equation in a maximally symmetric space 
except that the " Gravitational constant" is G(l — A/4 — r/6) _1 . When this 
constant is oo we can perhaps expect a de-Sitter solution even in the vacuum 
case i.e. when T = 0.Q 



3 Robertson- Walker Case 

For simplicity we consider only the Ricci squared term so setting r = 0. We 
also take a flat k = Robert son- Walker metric: 

ds 2 = -dt 2 + a 2 (dr 2 + r 2 {d0 2 + sin 2 6d<f) 2 )) (14) 

We use the conventions of Weinberg [24] such that the Ricci-scalar R is related 
to the Hubble parameter H as, 

R=-12H 2 -6H (15) 

With this metric the field equations give just two independent equations. 
The trace equation 

§-!)* + 24H> + 396^ + 2016§! + 9M* * + 936ff=| 

+216ff 4 |j - 72H'^ = 5^ (p - 3p) = k, (16) 
and the (0, 0) component equation 

-^P— (17) 

We have also included a perfect fluid matter source with equation of state 
p — (a — l)p (p and p are the pressure and energy density respectively). The 
conservation equation is given by 

p = SaHp ^P = ^t ( 18 ) 



4 As pointed out by the referee this argument is somewhat unsatisfactory, but does 
give a heuristic argument as to why de-Sitter solution are present. 
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A few well known equations of state are a = 4/3 (radiation); a — 1 (matter) 
and a = (cosmological constant) 

Scaling properties of equations (16) and (17) with respect to changes in 
the time scale suggest we look for solutions for which R and H 2 are propor- 
tional to each other i.e. to look for solutions for which 

a = -Ji (19) 

where a is a constant to be determined. The —ve sign is chosen to keep a+ve 
in our notation, a plus sign would require a large velocity H and would not 
be expected to be inflationary cf. eq.(15). Using the relation eq.(15) and the 
ansatz eq.(19) gives for the Hubble parameter 

H = 1 + H {2 - l/6*)t (20) 
Since a = aoexp(J Hdt) this leads to solutions of the form 

a(t) =ao(l + ^o(2-— )t) 12CT " (21) 

with H D and a Q the initial values of the Hubble parameter and scale factor 
respectively. Eq.(21) therefore includes solutions like a ~ t m with m = 

6(7/(12(7-1). 

If (2 - 1/6(7) = then simply H = H and 

a = a exp(tH ) (22) 

This usual de-Sitter solution corresponds as expected to a = 1/12 when 
R = -12H 2 . 

The ansatz when substituted into the Trace equation and the (0, 0) com- 
ponent equation , gives respectively 

R ( 12a 2 + Aa - - + \) = - Kl (23) 



Vacuum Case 



3 A 

3aR (l2a 2 + 4(7 - - + j ) -h 2 (21) 
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We first consider the absence of matter i.e. K\ — K2 — 0. Since we require 
R 7^ we set the a expression to be zero. 

12a 2 + 4a~l + \ = (25) 



(26) 



3 A 

with solution 

1^3/ V 

<7 = -6 + l-( 1 -A, 
The other solution of the quadratic equation can be ignored since we 
require a + ve. This gives for m 

m = ^ ^ (27) 

-3 + 2V3(l - \) 

A positive: 

Real i.e. non oscillatory solutions, require A > 1 ; the de Sitter solution 
m — > 00 occurs when 3 = 2-^/3(1 — 1/A) 1 / 2 =>- A = 4. We can plot the 
behaviour of m against A in Fig.(l) . For A > 4 we have power law inflationary 
behaviour with m ~ 1.6 as A — > 00. 

We consider the stability of these solutions with respect to linearized 
perturbations in the appendix. We find that they are stable to such pertur- 
bations. 

For 1 < A < 4 the universe is contracting i.e m < 0. There is also an os- 
cillation region for < A < 1 where m becomes complex: but since we are 
primarily interested in inflationary behaviour we do not consider it further. 
A negative: 

As A decreases the expansion rate m increases with m ~ 1.6 as A - >— 00 
. To get Inflationary behaviour (i.e. m > 1) requires A < —3. Again these 
inflationary solutions are stable to linearized perturbations. 
Matter Case 

We now consider the addition of an energy density such that 

P = ^ (28) 

We are setting n = 3a cf. eq.(18). Sticking with our same ansatz the trace 
equation and (0, 0) equations gives respectively, 

g-l-^-^Wgffi*' 4 - 3 "' (29) 
3 A / c 2 A a n v 1 
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2 1 te -12^« = ^A (30) 

3 A J c 2 A 3aa n v ; 

Equating the two gives 

h = A ~ 3a (31) 

The solution has the expansion behaviour m = 2/n = 2 /3a which is the 
usual relationship in standard General Relativity. 
Since for the solution eq.(21) we have 

(\ l/m 
m (32) 

the Ricci scalar can be written as 

R = --Hy /m a(t)- 2 ' m (33) 
a 

where the initial values H Q , a Q and p a are related by eq.(30). 
We have to avoid R = which occurs when \ja — i.e. 4 ^ 3a or n ^ 4. 
This is the equation of state for radiation whose scale factor behaves as 
a ~ t 1 / 2 . This is not allowed in this theory but a correction a ~ t 1 / 2 + 5 with 
S small is. 

Cosmological constant 
This is a special case of the previously example with n = i.e. p = p . Since 
we expect de-Sitter solutions we consider H = constant. From the Trace and 
(0,0) equations, we get 

12 (i - 1) " 2 = Kl (34) 

L 1 
A ~ 4 

Therefore 



I)" 2 'V CM) 



4 = 51= !_??!U P =-p (36) 

«2 V P / 

A de-Sitter solution occurs for the equation of state p = —p or equiva- 
lently a cosmological constant A ~ 8iiGp 

Because H should be real the value of A is constrained, from eq.(35) 

therefore A < 4 is required for such solutions. 



8 



4 Discussion and Conclusions 



We have found that provided the value of |A| is sufficiently large the alter- 
native higher order Lagrangian can have power law inflationary solutions, 
which are stable to linearized perturbations. Once matter is introduced, per- 
haps by a quantum process (c.f. particle production see e.g ref. [25]) in the 
early universe, the usual behaviour of General Relativity(GR)results. This is 
provided the matter is allowed to dominate while the tendency of the vacuum 
inflationary solution is to dilute any matter present. If instead quantum pro- 
cess proceed abruptly it might cause the inflationary phase to rapidly stop-a 
graceful exit problem. Until the coupling to matter was known it would be 
difficult to calculate the effects of particle creation and how such a scenario 
could proceed: whether such matter creation can cause the inflationary be- 
haviour to gradually switch over to a normal expansion a ~ t p with p < 1. 
In any case, there have been claims that the usual inflationary models are 
unstable to quantum processes [26]. If this is the case then not having a 
"classical" ending to inflation need not be a drawback: inflation would end 
due to quantum instabilities and a mechanism to enable a switch to a clas- 
sical non-inflationary phase is provided. A weak quantum instability would 
therefore appear to be a necessityin this, and similar power law inflationary 
models-weak in the sense that there is still a quantum probability of sufficient 
inflation occurring. 

These quantum processes would also be required to heat the universe since 
there is no oscillating field present: this is also a problem with scalar field 
power law inflationary models where the field potential is too shallow. If such 
processes are not present or cause the inflationary expansion to prematurely 
stop then a second period of inflation would be required. The first period of 
inflation would still have contributed to the flatness and horizon problem so 
that the second inflationary stage (presumably due to a matter field) could 
be shorter than the usual ~ 70 e-foldings. 

There is the problem that the radiation expansion a ~ t 1 ^ 2 cannot quite oc- 
cur (this would set the Ricci scalar zero). This might conflict with the usual 
predictions of nucleosynthesis and so require that A be smaller .| 

5 Or at least enabled the universe to wait, without recollapsing prematurely c.f. ref. [27], 
till a lower energy scale inflationary mechanism could dominate. 

6 We have not yet calculated how A affects the nucleosynthesis rates and if any possible 
differences disappear as A — > 0. 
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We could still envisage that during the early epoch of the universe the 
value of A is large: so giving inflation. As quantum gravity effects weaken due 
to the expansion of the universe the value of A reduces, so effectively agreeing 
with GR later. We would hope to show that as A — > we would recover or 
agree with the results of GR. In this case A varies like a "running coupling 
constant" in particle physics due to the relevant energy scale present. 

An alternative scenario would be for A to be simulated by the effect of 
extremely highly energetic matter e.g. a string epoch. A could then change 
suddenly rather like the induced gravity model c.f ref.[23] due to the dif- 
ferent behaviour of the matter fields present (in the string epoch this could 
correspond to a phase transition ). This would require us to understand the 
matter fields present, and their coupling to the higher derivative terms, be- 
fore we would know if they could set A small. Note that the field equations 
have been derived under the assumption that A is independent of time. If this 
assumption was dropped then the field equations would be more complicated. 

Because the effective Newton's constant G e ff can change sign there is also 
the chance of finding wormhole solutions. This is similar to the wormhole 
found in the case of a conformally coupled scalar field which has G e ff < 
in its vicinity [28, 29]. Wormholes can be found in the usual higher derivative 
Lagrangian (1) if we take the "wrong sign" of epsilon[30,31]: this is related to 
the bounce (avoidance of the singularity) in such models for spatial curvature 
k = -1[31]. 

In order to find wormholes (or bounces) in the Lagrangian presented here we 
would have to keep the spatial curvature k in the equations. This would be 
rather more complicated as would the case of Bianchi models: these could 
perhaps be more realistically tackled using algebraic computing packages. 

In conclusion we have obtained vacuum inflationary solutions for a non 
analytic higher order Lagrangian. This Lagrangian might contain certain 
aspects of the "correct quantum gravity" theory at high energy densities. 
Such solutions are further evidence that De Sitter or power law inflationary 
solutions are a general feature of higher order gravitational theories. 

Acknowledgement C.Xu would like to thank Prof. P. Havas for his 
enlightening discussions on various aspects of the alternate gravity theory. 
We are grateful to the referees for many useful comments and suggestions 
which are reflected in this revised version. 
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Appendix 

The question of stability of the solutions is clearly of importance. However, 
a comprehensive discussion of this topic would include such concepts as e.g. 
linearized, weak, asymptotic, and Liapunov stability. This would entail con- 
siderable work beyond the scope of this paper, so we limit ourselves to just 
an analysis of the linearized stability of the vacuum equations. For a general 
reference see eg. [32]. 

First consider the (0,0) eq.(17) rewritten in the form 

together with the expression for the Ricci scalar 

H = -2H 2 - R/Q = F 1 (H, R) (39) 

We consider perturbations h and r around the known solutions such that 
now 

H = H + h R = R + r (40) 

where H and R are the previously found solutions. The linearized equations 
for the perturbations are given by 

h ^j^ h + aj ^ r (41) 
f = aj 2 ijm h + ej^B lr (42) 

Calculating the partial derivatives eventually leads to 

h = -AhH -r/6 (43) 
f = fchH 2 + (3 2 rH (44) 

H' 1 = Hq 1 +yt, 7 = (2 - l/6<r) > (45) 



where 
and 



A = (rTW( 18+lte+ i + ^.) (46) 
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02 = 



(1 + 3(7)2 V 2 2 



1 45 



2 \ 



-~ - — (x-45(7 z + — • 



3(7 y 



(47) 



In deriving /3i, /3 2 we have made use of eq.(25) to express A in terms of 
(7. Scaling to a new time parameter r such that dr/dt = 7 and renaming 
u(t) = h(t) and v{r) = r(t) gives 



. . 1 / 4u v 
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1 (fan I3 2 v s 



(48) 



(49) 



We have to find solutions with r(= H 1 ) > r = H 1 , with initial conditions 
u(t ) = m and d(t ) = t>o- 

The nature of the equations suggests we look for solutions of the form 



U = UqT a V = VqT a 1 



(50) 



When substituted into the equations this ansatz leads to the linear system 
for u ,v 

(4 - a-f)u + v /6 = (51) 
Piuo + (02 + (a + 1)7) ^o = (52) 
Non-trivial solutions require that the determinant vanish i.e. 



or 



a 



(4 - cry) (ft + (a + 1) 7 ) - 0i/6 = 



with solutions 



«1,2 



4-0 2 

7 



7 



- 1 



± 



\4 



4-02 

7 



- 1 



_ 40 2 -A/6 4 

7 2 7 



(53) 
(54) 

(55) 



which can be written in the form (with obvious definitions for A, B) 



C*l,2 = - ± -VA 2 + 5. 

7 7 



(56) 
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For stability we require both oi and a>2 > this occurs if the following 
inequalities are satisfied (recall 7 > 0): 

A>0,B<0,A 2 + B>0 (57) 

Substituting for fiijfa allows these inequalities to be written as 

< (TtW (I + i6a + 63,72 ~h) (58) 
0< (rrW( 1 + 51ff+108ff2+ i^-^) (59) 

1 / 14 2 \ 

0<(4 -*^» + (TTw(" 4 " 2 ° 4 ^ 432ff + S-9^) ' (60) 

The value of a when inflationary behaviour occurs i.e. A > 4 is constrained 
in the region 

i<«,<QV3-i)=. 0.122 (61) 

where we have again used eq.(26) to relate a to A. For A < — 3 the corre- 
sponding expression is 

±(V3-1)<<7<± (62) 

With these values of a the inequalities are satisfied. The first two are eas- 
ily satisfied while the third one has been checked numerically. In fact, the 
restriction to the inflationary regime is not strictly necessary since the in- 
equalities are only violated in the oscillatory region < A < 1. 
The values of 01, a 2 are therefore positive and the perturbations decay as t or 
t — > 00. Let us represent the solutions for the perturbations as cu(ai), 07(0:2) 

, where to = ^ ^ ^ . These two independent solution of u will span the so- 
lution space of eq. (48,49), which is a two dimensional vector space. This is 

H 
R 



in contrast with the solutions (which can be called Q = ( ^ ) ) which go as 



H = 1/r, R = -1/(ctt 2 ) (see eq.(19). 

In order for the perturbations not to dominate they must decay more rapidly 
than their corresponding solutions i.e. we require a > 1. Because o 2 is less 
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than Oi, 00(0:2) will give the dominant behaviour as r — > 00. 
Now «2 is given by 

o 2 = - - -VA^TB. (63) 

7 7 

Numerically it can be shown that (A 2 + B) ~ 0.2 so approximately 

o 2 ~ - = -— (- + 36a + 63a 2 - — ) - — (64) 

7 2 7 (l + 3a) 2 V2 2a/ 7 K J 

For a in the inflationary range (A > 4) we obtain 

~ 1.4 < o 2 < 00 (65) 

(For A < — 3 the corresponding expression is 1.4 < a 2 < 1.8.) 

We can therefore conclude that since a 2 > 1 the perturbations 07(0:2) will 

decay faster than their corresponding solutions Q. 

Again this restriction to the inflationary range was not strictly necessary 
except for the oscillatory region. Although we do not try to analyze this 
region further, we would also hope that for stability in this region where a 
is complex, the real part i.e both ot^ ai are positive. 

We would also like to show at some stage that the inclusion of matter is 
likewise stable. Although we are fairly confident that this will be the case, 
since more rapidly expanding solutions are usually more prone to instability, 
it would require extensive further work to be certain. 

The trace equation should also be considered but it can be shown that it only 
imposes restrictions on the initial values Wo^o an d so does not add further 
to the stability analysis. We can therefore conclude that the inflationary 
solutions are stable, at least to linearized perturbations. We can be fairly 
confident that other solutions (i.e oscillatory, matter dominated) will also be 
stable without explicitly doing each case. 
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Figures 

Fig.l 

Plot of the expansion coefficient m against the parameter A. Power law 
inflation occurs for A > 4 and A < — 3. 
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